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It has recently been found that supernova explosions can be simulated in the labo- 
ratory by implosions induced in a plasma by intense lasers. A theoretical explanation is 

>' 

Q\ ' that the inversion transformation, (E : t — 1/t, x — > x/i), leaves the Euler equations 
ON ! 

^ | of fluid dynamics, with standard poly tropic exponent, invariant. This implies that the 
O 

kinematical invariance group of the Euler equations is larger than the Galilei group. 
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In this paper we determine, in a systematic manner, the maximal invariance group Q 



of general fluid dynamics and show that it is a semi-direct product Q — SL(2, R) A G, 
where the SL(2, R) group contains the time-translations, dilations and the inversion E, 



and G is the static (nine-parameter) Galilei group. A subtle aspect of the inclusion of 
viscosity fields is discussed and it is shown that the Navier-Stokes assumption of con- 
stant viscosity breaks the SL(2,R) group to a two-parameter group of time translations 
and dilations in a tensorial way. The 12-parameter group Q is also known to be the 
maximal invariance group of the free Schrodinger equation. It originates in the free 
Hamilton-Jacobi equation which is central to both fluid dynamics and the Schrodinger 
equation. 
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Introduction 



Considerable efforts are being devoted at present in centres such as the National 
Ignition Facility at the Lawrence Livermore Laboratory, U.S.A., and the Laser Mega- 
Joule in Bordeaux, France, to simulate astrophysical systems in the laboratory. The 
motivation for this programme comes from observational evidence [1] that the dynamics 
of the mixing of gases during a supernova explosion [2] is very similar to the turbulent 
splashing and mixing of a plasma as a fusion capsule is bombarded by high intensity 
laser beams [3]. Since the former physical situation deals with an explosion while the 
latter concerns an implosion, and because the time and length scales involved in the 
two cases are drastically different, it was a challenging task to produce a theoretical 
explanation for the observed similarity. 

In a recent paper [4] , Drury and Mendonga produced such an explanation by show- 
ing that the Euler equations of fluid dynamics which govern both the systems are left 
invariant by an inversion transformation E : t — > — 1/t, x — > x/t. The discovery of 
this transformation was inspired by similar transformations of the form x — > /(£)x that 
are used by cosmologists to factor out an arbitrary uniform expansion or contraction of 
the system [5]. As is well-known, the Euler equations of fluid dynamics are invariant 
under the Galilei group [6]. The results of [4] show that the maximal kinematical 1 in- 
variance group of fluid dynamics Q is larger than the Galilei group [7] . In this paper we 
determine, in a systematic manner, the structure of Q for general fluid dynamics. We 
also discuss subtleties associated with the inclusion of viscosity fields and show that the 
Navier-Stokes assumption of constant viscosity leads to a reduction of Q in a tensorial 
manner. 

1 Meaning transformations involving only space and time coordinates as opposed to 
those involving internal or gauge degrees of freedom. 
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It will be shown that Q is of the form 

Q = SL(2, R) A G (1) 
where G is the connected, static Galilei group which induces the transformations 

x -> Rx + a + vt, t->t (2) 



and SL(2, R) is the group 

at + (3 



aS-p^ = l (3) 



which includes parity, time translations, dilations and the inversion E. 

Q is isomorphic to the Niederer group, which is the maximal invariance group of 
the free Schrodinger equation [8]. The reason for this isomorphism is that the kinetic 
part of the fluid-dynamic Lagrangian (see (17)) is of the free Hamilton- Jacobi form, 
which is the classical limit of the Schrodinger operator. 

The Strategy of Proof 

The general fluid dynamic equations in n-dimensional space are [9] 

Dp = -pV • u (4a) 

pDu = -Vp + V (46) 
De = -(e + p)V • u (4c) 
where the convective derivative D and the viscosity terms V are defined by 

D = — + u • V and Vi = Vj (rj(V jUi + V l u j - -5ijV k Uk) ) + V^CVfc-Ufc) (5) 

respectively. In the above equations p, u, p and e stand for the density, the velocity 
vector field, the pressure, and the energy density of the fluid respectively and 77, ( are 
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the shear and bulk viscosity fields. The above differential equations are augmented by 
an algebraic condition called the polytropic equation of state which relates the pressure 
to the energy density thus: 

V = (lo - l)e =^ V + e = 7oe (6) 

7o being a constant called the polytropic exponent. This equation may be used to 
eliminate p from (4). Further, by making the substitution, 

e = XP 10 (7) 

the original set of equations (4) reduces to 

Dp = -pV • u (8a) 

pDu = _( 7o -l)V(xp 7o )+V (86) 
D X = (8c) 

To find the maximal invariance group of these equations we first note that, for any sys- 
tem (such as the above) for which there is no feedback from the fields to the space-time 
transformations, the maximal invariance group Q of the general class of configurations is 
also an invariance group for any special sub-class of configurations; though not necessar- 
ily a maximal one, since a restricted set of configurations could have a larger symmetry 
group. Formally, 

c s cc^gcg s (9) 

where C denotes general configurations and C s denotes a special class. Thus a general 
strategy for finding maximal invariance groups is to first find a (tractable) sub-class of 
configurations C s for which the maximal invariance group Q s can be found, and then see 
what conditions are imposed on Q s by the general configurations. This is the strategy 
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we shall adopt; except that in our case we choose a sub-class for which the generalisation 
to arbitrary configurations imposes no further conditions i.e. Q s = Q. The tractable 
sub-class we choose is obtained by making three simplifications, namely, ignoring the 
effects of viscosity by tuning the viscosity fields r\ and ( to zero, by setting % = 1, and 
by letting the velocity vector u be curl-free. The resulting sub-class is the configuration 
of inviscid, isentropic and irrotational flows which are well-known in fluid dynamics. 

The Sub-class of Inviscid, Isentropic, and Irrotational Fluids 

The first simplification we consider is to ignore viscosity effects in which case (8) 
reduces to 

Dp = -pV • u (10a) 

P Du=-{ lo -l)V{ X p 1 °) (106) 

D X = (10c) 

Second, we note that it is consistent with the field equations to set % = 1, with u not 
necessarily curl-free, in which case the equations in (10) reduce to, 

Dp = -pV • u and Du = -7 (7 - l)p 7o " 2 Vp (11) 

The reason for making this simplification is that the system becomes a local Lagrangian 
one. This can be seen by using the standard Clebsch parametrization [10] 

U= _V0-*A70 (12) 

for the vector-field u, in which case the local Lagrangian density is, 

C = pty + vO - ^(V^ + z/V6>) 2 - p 70 " 1 ] (13) 



It will be convenient to write this Lagrangian density in the form 

dt 



C = p U--(V(l)) 2 + uV6 -— (V6) 2 - p^' 1 ] where V = ^ - V0 • V (14) 



is the convective derivative D restricted to the curl-free part of u. Clearly v is a Lagrange 
multiplier, and, by varying it, we obtain 



V9 



and L = p[j> - i(V0) 2 - \^ - p^ 1 } (15) 



(W) 2 ~ ~ ~ T 2 V ^ 2(V0) : 
Note that the first equation here is equivalent to D9 = since 



bq = vq- (uve) -ve = ve- u{vef = {vef - u ) ( 16 ) 

We now proceed to the third and last simplification mentioned. 

The third simplification we make is to consider the curl- free case u = — V0 or 
v = 9 = 0. Then the Action corresponding to the Lagrangian density (13) reduces to 

S = J d n xdt p - ^(V</>) 2 - p 7 ° (17) 

Note that the quantity in square brackets is just the Hamilton-Jacobi function for a 
free particle. Following our strategy we now seek the maximal invariance group of the 
Action (17). 

The most general transformation involving the fields, as shown in the Appendix, 
takes the linear, inhomogeneous form 

& = &(M)» T = r(x,t), <j> = s(x,t)<f> + \(x,t), p = fj,(x,t)p + fj,'(x,t) (18) 

where s is a constant. However, one sees at once that if p! ^ the form invariance of 
the Action is violated by inhomogeneous terms. Furthermore, since 

d dr d d^i d d dr d d^j d 

dt dt dr dt <9£j dxi dxi dr dxi d£j 



we see that we obtain a term (dcft/dr) 2 unless 

— - l 

dxi 
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(20) 



This means that the d<j)/dt term does not pick up an x- dependent factor relative to 
d<p/dT. Finally we note that under the transformations (18), the second term in the 
Hamilton- Jacobi part of the Action (17) gives 

1 d<t> d(f) 1 dcf) d(f) <9& d£j _ ( . , . , . . 

2 h -d^d^ ~* ^WiWj ere 9ij = dT k d^l hkl = Q{x)9i ^ (21) 

where is the flat Euclidean metric. As is well-known, the most general transfor- 
mations that have the above property are the special conformal transformations with 
parameter 6, for which O = (1 + 2b ■ x + b 2 x 2 ) 2 , and the Euclidean transformations for 
which = 1. It is easy to see however that, for the Action in (17) to remain invariant, 
the d(f)/dt part of the Hamilton- Jacobi function must transform in the same way as 
the (d(f)/dx) 2 part. But since (20) forbids d<j)/dt to pick up an ^-dependent factor, the 
conformal transformations in (21) are not allowed. Thus the possible transformations 
reduce to 

& = Mijtyxj + wtit), T = r(t), <j> = s<f> + \(x,t) p = n(x,t)p (22) 

where s is a constant and M is Euclidean i.e. M = f(t)R(t), where / is a scale factor 
and R(t) is a rotation matrix. In that case 

ii M tj ) ^ (detM^o) (23) 



/ <9r 


dr_ \ 




/ dt 


dxj 








H 


\ dt 


dxj ) 





and the Action becomes 



d n £dT 



detM^r 



. d ■ d (0 - A) detM 2 80 d\ 



f d n £dT r d tiWi^fl u detM 2 dj> d\ 2 g^" 1 )^" 1 
Invariance then requires that 

\2 _ Q + QT1 J /, -1\7 -1 



(24) 



sp(detM) 71 = 1, (detM) 2 = sf, and (/U -1 ) 70 " 1 = - (25) 



and 

( ^ + a7^ )((/, - A) -2br^ = a7~2 ( W } (26) 

Equations (25) can be readily solved to get 

2 

s = 1 =>■ (detM) 2 = f , // = (detM) _Tl , and 7o = 1 + - (27) 

n 

The result for 7 G is in agreement with [4]; in particular in three-dimensional space, 
n = 3, and 7 G = |, as expected. Substituting these results in (26), and requiring that 
it be satisfied for all d(j)/d^i, splits it into 

<9& dA , OX l/dA\ 2 2 

Notice that the second equation is just the free Hamilton-Jacobi equation with Action 
A. The maximal invariance group can now be obtained by solving the above equations 
for the functions £j and A. Thus the invariance group Q is determined by the structure 
of the Hamilton-Jacobi equation for a free particle. However, we will now show, even 
before solving these equations, that the group of invariance remains the same when the 
simplifications made in this section are relaxed and we return to the case of the general 
fluid configurations. 

Proof of Maximality for General Fluid Configurations 

We first note that the first and last equations in (28) are exactly the conditions for 
the restricted derivative V to transform covariantly since 

T>( x t 6) = — - ^tJL = (fJL + d& 9 \ _ ( d tj d£fc \ 94> g 
dt dxi dxi V dr dt <9£j / V dxi dxi J d£j d£k 

- (4 - -41) + (t — 2 |)|- - «^^.» 

(29) 



Let us now relax the assumption that u is curl-free, with % = 1 still being true. From 
the Lagrangian density for this case, 

£ = , [4 _l(V,)»-igg (30) 

and the covariant transformations, 

V(x,t,c/>) = (detM) 2 £>(£,r,0) and V(x) = (detM)V(O (31) 

we see that (30) will be invariant under the given space-time transformations provided 
only that v and 9 are scalars. Thus the case of general u puts no restrictions on the 
space-time transformations but determines the transformation properties of the v and 9 
components of u. In particular we see that the full convective derivative D and u have 
the transformation properties 

D(x,t,cf>) = (detM) 2 D(^r^), and u = _-L_ (u - VA) (32) 

Note that the inhomogeneous part of the u-transformation comes from the curl-free part 
of u. Furthermore we see that the Q transformations do not preserve the condition of 
incompressibility (V • u = 0) and hence we need to consider general, compressible fluid 
configurations. 

All this is within the Lagrangian framework. In order to relax the % = 1 simplifica- 
tion, we need to allow for general fields x- F° r this we must go outside the Lagrangian 
framework and consider the field equations. But that is simple since the field equations 
with % = 1 are invariant and we see by inspection that the equations for general x 
remain invariant provided that x transforms as a scalar. Thus the set of equations in 
(10) is invariant with respect to the maximal invariance group provided u, 9 and x 
transform as scalars. Thus the maximal invariance group for v = 9 = 0, x = 1 (cor- 
responding to our sub-class without viscosity) is the maximal invariance group for all 
inviscid configurations. 
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Finally we include viscosity effects. It is straightforward to verify that the viscosity 
terms in (4b) are covariant only if V • u = or V^C = 0. However, as already mentioned, 
the former condition is not invariant under Q transformations and hence covariance 
requires that ( is a constant in space (although not in time). In this case, the fluid 
equations (4) would remain invariant if the viscosity fields transform as tensors of rank 
n, where n is the space dimension i.e. r\ — > detM n ?y and ( — > detM n (. This in turn 
implies that the equations would remain invariant under the Navier-Stokes assumption 
of constant viscosity only if detM is a constant which breaks the Q symmetry in a 
tensorial way. Therefore the maximal invariance group of the Navier-Stokes equations 
contains just the static Galilei group plus dilations and time translations corresponding 
to detM being a constant, as will be clear when we solve for the Q transformation 
functions in the next section. 

It is worth mentioning that the breakdown of Q by the Navier-Stokes restriction 
does not contradict (9) because, in contrast to the other restrictions 1] = ( = 0, x = 1, 
V x u = that we have considered, this restriction has a feed-back effect on the non- 
linear part of the space-time transformations. 



Returning to the equations in (28), it is first useful to note from the definition of £ 
in (22) that 

= Aidj+Wt where A = M T M T = f T f~ l +R T R T and W = M^-(M _1 «;) (33) 
or or 

where the subscript r denotes differentiation with respect to r. Hence the first equation 
in (28) becomes 



The Solutions for the Transformation Functions 




(34) 
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Here the equation in the brackets shows that A is a symmetric matrix. But since R 
is a rotation matrix, R T R T is in the Lie algebra of the rotation group and hence is 
antisymmetric. Thus we have the result R T = and hence R is a constant (rigid) 
rotation matrix. 2 

It then follows that detM = / and from (28) and (32) we have 

t = f and u = - (u - VA) (35) 
The definitions of A and W then simplify to 

A = f T f- 1 and Wi = fd T (^j) (36) 
Integrating (34) with respect to £ we have 

A = -A^ - W& - h(r) (37) 

Substituting this into the second equation in (28) we get 

(0 2 dA BWi dh \, AJ . Tir . 2 
-f" = 2^ W) (38) 

Comparing the coefficients of the powers of £ breaks this into 

OA a2 dW l . <9 T f d/i W 2 

8F = " A - -g? = -W = -f-n. Tt =~— (») 

We can solve the second equation in (39) explicitly to get 

Wi = J d t (^pj = Vi => Wi = f(t)( Vi t + a t ) (40) 

where the i>; and are constants. It follows that 

i% = fxi + Wi = f(t)(xi + a i + v l t) (41) 



With this in mind, we suppress the matrix for convenience, and restore it later. 
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The third equation in (39) can then be written as 

dh v 2 dh v 2 , , v 2 

Since the transformations of the fields are given by p = f~ n p and = + A, we see 
that everything is determined by the function / which satisfies 



d ff\ ( f 



dr\pj V/3 
Using f = / 2 , this is easily seen to be equivalent to 



2 



(43) 



d 2 f = and L-l(Z.) 2 = o (44) 



f 2 Vf , 

The left hand side of the second equation is the Schwarzian derivative of r and it follows 
that the general solution of the above equation is 

t = — ~^-^r and f = where a5 — #7 = 1 (45) 

Thus the most general transformations of the time coordinate form an SL(2, R) group. 
The second equation in (35) then reduces to 

d = ( 7 t + 5) u _ 7 (x + a) +5v (46) 

From (41) it also follows that 

& = ^^-^(x, + a, + v.t) (47) 

Returning to the general case, and using the inversions 

— 1 , 7 (a; + a) — Sv . A . 

1T - C ' = —S md TtTi (48) 



we have 

(7^ - 1>) 2 , [7(2; + a) - <k>l 2 , , /, 5v 2 \ 



27(77" — a) 27(7^ + 5) \ 7 
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The case detM = / = a constant corresponds to letting 7 = and from (45) and (47) it 
is clear that the group of transformations in this case includes the static Galilei group G, 
the dilations, the time translations and parity, but excludes time-reversal. As already 
mentioned, this is the maximal invariance group of the Navier-Stokes equations. It is 
also useful to consider the following two special cases. 

I. j3 = 7 = 0, a = 1: The Connected, Static Galilei Transformations: In this case, we 
have 

g : r = t, £ = Rx + a + vt (50) 

where we have restored the rotations, and from (49) it follows 

v 2 

A = -h - v • (Rx + a) 1 (51) 

These equations describe connected, static Galilei transformations which exclude parity 
and time-reversal. 

77. a = v = 0, R = 1: The Inversion Transformations: In this case, we have 

"■■ aS -^= i < 52 > 

-yt + 7t + o 

These are the SL{2, R) generalisations of the inversion transformations presented in [4]. 
For the transformations of the fields we have from (49) 

x= 2WTs)- ho (53) 

The Maximal Invariance Group Q 

To understand the structure of the group, we study the relationship between the 
SL(2, R) group and the connected static Galilei group G. Let us first consider a conju- 
gation ofajeGbyaae SX(2, R). By making three successive transformations of x 
and t we find that 

a' 1 (a, p, i)g{R, a, v)a(a, (3, 7) = g(R, a a , v CT ) (54) 
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where 

5 0\ fa 



Vcr j v 7 a/yv 
This shows that G is an invariant sub-group and so the group structure is 



(55) 



Q = SL(2,R) AG (56) 

where A denotes semi-direct product with G as invariant subgroup. More precisely, if 
we recall that G itself takes the form 

G = R A (T(a) (8) B(v)) (57) 

where T and £? are the translation and boost groups with parameters a and v respec- 
tively, then we see that SL{2, R) commutes with R and mixes T and B in the manner 
shown in (55). 

The original inversion E is the special element of SL(2, R) for which (a, /?, 7, 5) = 
(0, —1, 1, 0). Note that T? = P where P is the parity. Furthermore if we consider the 
coset elements g ^ (R, a, v) = T,g(R, a, v), where g G G, we have using (54) and the 
standard property of Galilei transformations viz. g (R, a, v)g(R', a', v') = g(RR' , Ra' + 
a, .Rv' + v) , 

gjR!, a', v')g s (R, a, v) = g p (R'R, R'a-V, R>v + a') 



(58) 

=> g\ (R, a, v) = g p (R 2 , Ra - v, Rv + a) 

where we have used the obvious notation g (R, a, v) = Pg(R, a, v). It follows from 
the above equation that 

g\{R, a, v)=g(R\ (R 2 - l)(Ra - v), (R 2 - l)(Rv + a)) (59) 

Since Ra — v and Rv + a are linearly independent, this shows that every connected 
Galilei transformation is the fourth power of a coset transformation. 
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Extension to the Quantum Theory 

If a classical system is described by the Galilei group G in (57) then, in the field 
theoretic representation, there is a central extension of G; namely, the well-known one- 
parameter mass group whose generator commutes with all the generators of G. Since 
a quantum wavefunction can be thought of as a nonrelativistic field, the corresponding 
quantum system is described by a central extension of the group of the form 

G = RA (r(a),B(v)) (60) 

where the group (T, B) is no longer abelian, but is a Heisenberg-Weyl group of the form 

T(a)S(v) = B( V )T(a.)e lhMa - v (61) 

and M is a central constant. It is obvious that this relation is invariant with respect 
to the transformation (55) of the parameters a and v induced by SL(2,R). Thus Q 
remains a symmetry group of the quantised system. 

Connection to the Free Schrodinger Equation 

The fact that Q is a covariance group of the Hamilton- Jacobi function d<j>/dt + 
l/2(V</>) 2 shows that it is an invariance group of a free particle. The above discussion 
of the quantum extension implies that it is also an invariance group of a quantised free 
particle. Indeed it was shown in 1972 that Q is the maximal invariance group of the 
free particle Schrodinger equation [8]. Since the invariance under the Galilei group is 
well-known, it suffices to verify that the Schrodinger equation 

cty n 2 d 2 ^ 

th dt + 2^d^ = ° (62) 
remains form-invariant under the Q. It is easily checked that this is accomplished by 
the following transformation of the wavefunction 

^(x,t) ~ (7T-a)^e- iA V>(£,T) ~ (~ft + 5)-%e- iX if; (63) 

15 



where A is defined in (49). Note that at time r = oo, i.e. t = — ^, the wavefunction 
becomes infinite, but this is precisely the singularity associated with the explosion. 

Conclusions 

We have investigated the maximal kinematical invariance group of general fluid 
mechanics and have found that it is a semi-direct product of the form SL(2,R) A G 
where G is the static Galilei group. The incompressibility condition is not preserved 
by the above transformations and hence the viscosity terms transform covariantly only 
if the viscosity field Q is constant in space (although not in time). In this case the 
fluid dynamic equations remain invariant if the viscosity fields transform like tensors 
of rank n, where n is the space dimension. The Navier-Stokes assumption of constant 
viscosity breaks the SL(2,R) part of Q to a two-parameter group of dilations and time 
translations in a tensorial manner. 

The inversion transformation E found in [4] is a special element of the SL(2, R) 
and acts like the square root of parity. The transformations generated by the coset 
elements g = T,g act like the fourth roots of connected Galilei transformations. It is 
also pointed out that Q is the Niederer group which is the maximal invariance group of 
the free Schrodinger equation. The reason for this is that Q is actually the invariance 
group of the free Hamilton-Jacobi function, which is common to both fluid dynamics 
and the classical limit of the Schrodinger equation. It is surprising that such a basic 
result does not seem to be more widely known! 

Appendix 

In this appendix we show that the linear inhomogeneous transformations of the field 
variables considered in (18) are indeed the most general field transformations allowed. 
In order to show this we begin by noting that the most general transformations can be 
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written as follows: 

4> = 4>(£,t,<I>) and p = p(£,r,p) (Al) 

where 4> and p are a priori arbitrary functions of the old fields 4> and p. Note that we 
are not allowing the fields to mix, but this is reasonable since 4> and p are really not on 
the same footing; whereas <fi appears in the Action (17) only through its derivatives, p 
appears without any derivatives. Using the fact that x, t — > £, r, the Hamilton- Jacobi 
function can be rewritten as follows: 

dt 2 K dx> [ dtd£ dtdr> 2 [ dxd£ dx dr' [ ' 

Now it follows from 

4> = 0(C,r,0) cf> = F(£,T,$) (A3) 

that 

dcj) dF 04>dF d(f> dF 84>dF , . t . 

di di d£ d(p or dr dr d<p 

We would now like to fix the form of F by requiring the covariance of the Hamilton- 

Jacobi function i.e. 

d(j) 1/9^x2 d(j) 1/<9(/k2 

~dt~2^dx' K dr~2^W 

Substituting (A4) on the right hand side of (A2), we notice that we get terms of the 

form (d4>/dr) 2 unless dr/dx = =^ r = r(t). Equating the coefficients of d(f)/dr and 

— (l/2)(d(f)/d£ t ) 2 because of covariance, we get 

5 = # {M) 

Requiring the coefficient of the <9(/>/<9£ term to vanish produces 

6 i = # (A6) 
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Similarly requiring the coefficient of the (^-independent terms to vanish we get 

d£dF_ dr_dF_ 

Substituting this equation in (A6) gives, after some algebra, 

of__iAM_ (A7) 

\dx) at 

From (A5) it follows that since <p is arbitrary, the left hand side changes whereas the 
right hand side, for a given £, r, x, t, is a constant in the field. Hence 

where s _1 is a constant in the field variable. Solving the above differential equation we 
get 

F = s- 1 ^ + A(£, r) = r)0 + A(£, r) (A9) 

where A = — sX. Further by taking a derivative of (A9) with respect to £ we see that 

OF ds' 1 ~ ,84> OX . A . 

Again from (A6) it is clear that the left hand side of the above equation is independent 
of 4> and hence, s _1 is independent of £. Similarly, by appealing to (A7) we can show 
that it is independent of r. Hence s is a constant. A similar analysis may be carried 
out to establish the generality of the p transformations in (18). 
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